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'nI- , Abstract 
O 
O 

' Recent calculations of s-wave pion production have severely underestimated 

0^ ' the accurately known pp ppir^ total cross section near threshold. In these 
calculations, only the single-nucleon axial-charge operator is considered. We 

^ ' have calculated, in addition to the one-body term, the two-body contributions 



to this reaction that arise from the exchange of mesons. We find that the 



^ ' inclusion of the scalar cj-meson exchange current (and lesser contributions 



from other mesons) increases the cross section by about a factor of five, and 
leads to excellent agreement with the data. The results are neither very 
sensitive to changes in the distorting potential that generates the A^A^ wave 
function, nor to different choices for the meson-nucleon form factors. We 
argue that pp — > ppn^ data provide direct experimental evidence for meson- 
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exchange contributions to the axial current. 
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I. INTRODUCTION 



Recently, a precise measurement of the total cross section for pion production in the 
reaction pp —* ppix^ was carried out using the electron-cooled stored beam of the Indiana 
University Cyclotron Facility (lUCF) Cooler with an internal gas target [jl|. The novel 
technology made it possible to extend the measurements to within a few MeV of threshold. 
Over most of the covered energy range, contributions from higher partial waves are negligible, 
and the cross section is thus due to a single partial wave (^Pq ^ ^•S'o + s-wave pion). That 
this is the case can be deduced from the energy dependence of the cross section and the 
angular distribution of the outgoing protons 

As we will discuss below, s-wave pion production is sensitive to the axial charge of the 
two-nucleon system. A number of calculations using the single-nucleon axial-charge operator 
0J^,^,0] have been carried out with various choices for the A^A^ distorting potentials. Quite 
surprisingly, such calculations all underestimate the data by approximately a factor of five. 
This discrepancy is even more serious since there is little ambiguity in the calculations, 
because only a single partial wave is involved, and the A^A^ wave function and the (free) 
operator are both well known. 

The discrepancy could be explained by an enhancement of the axial charge in the A^A^ 
system, which could come from a relativistic effect. There exist many relativistic calcula- 
tions for nuclear systems; for example, relativistic impulse-approximation calculations [^,|^ 
accurately reproduce elastic proton-nucleus scattering data. Relativistic models character- 
istically feature large Lorentz scalar and vector potentials 0. The strong scalar potential 
enhances the lower component of the Dirac wave function of the nucleon; this change en- 
hances the axial charge |^ (see Sec. Several authors have examined the consequences 
of this relativistic effect on (3 decay or muon capture rates in A = 16 nuclei and nuclei 
near A = 208 . 

In nonrelativistic models, relativistic effects can be incorporated formally via meson- 
exchange currents (MEC's). In this view the relativistic effect would be represented by a 
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MEC involving a scalar a meson. Note that a exchange also provides a phenomenological 
model for the intermediate-range attraction in the A^A^ interaction and is responsible for 
part of the spin-orbit potential. Many authors have examined MEC contributions to the 
axial-charge operator in nuclei, and found that the largest contribution is due to vrp exchange 
TT|, which enhances axial-charge matrix elements by about 60%. The next largest MEC 



is believed to be the relativistic a exchange, contributing another 40% to the axial charge. 
Together these contributions could account for the approximately 100% enhancement seen 
in axial-charge matrix elements for a variety of j3 decays ||lT|. However, the size of the 
MEC contribution that one would derive from experiment depends on assumptions about 
complicated wave functions and short-range NN correlations. 

Because of isospin considerations, the vrp current (proportional to the dot product of the 
isospin matrices fi x r2 and the pion field) cannot contribute to pp ppii^ . However, other 
MEC may well be important. In this case, pion production in the pp system would provide a 
unique laboratory for testing MEC models for a number of reasons. First, the initial proton 
must decelerate in order to produce a pion. The resulting large momentum mismatch can 
best be mediated by meson exchange. Heavy mesons are favored; indeed, the initial relative 
momentum near threshold of about 1.9 fm~^ is comparable to the mass of the a meson. 
Thus the process is very sensitive to two-body contributions. Second, contributions from 
intermediate states with a nucleon and a spin-| A are expected to be small since such a 
system cannot be formed in a relative s state because of angular momentum considerations. 
Furthermore, the A cannot decay into a spin-| nucleon and s-wave pion. (Reference 0] 
provides estimates of A contributions.) Third, pion rescattering effects are small because 
the scattering length for pion-nucleon s-wave scattering in the required isospin combination 
is small, as will be discussed in Sec. |1[ Finally, the wave functions in this case are simple 
and can be calculated reliably. 

Lee and Riska ||12| have recently suggested that the inclusion of MEC could explain the 
magnitude of the observed pp ppix^ total cross section. They calculate MEC by assuming 
a simple operator form for the A^A^ potential. This allows the calculation of MEC for 



phenomenological potentials, but depends on model assumptions. In this paper, we perform 
a similar calculation; however, we use an explicit one-boson-exchange model for the A^A^ 
interaction as well as for the calculation of meson-exchange contributions. We also include 
the Coulomb interaction, which is quite important near threshold. Furthermore, we examine 
the sensitivity of the calculated cross section to many of the model ingredients. 

In Sec. H, we outline the formalism for the calculation and discuss the input parameters 



and some computational details. In Sec. |T| we collect and discuss the results for a variety 
of one-body and two-body contributions and for different NN interactions. Section lists 
the conclusions from our work. 

II. FORMALISM AND DETAILS OF THE CALCULATION 

In this section we present the formalism used to calculate the total pp — * ppn^ cross 
section for s-wave pion production. The calculation is carried out in coordinate space because 
this allows for a simple treatment of the Coulomb interaction between the two protons. The 
first three subsections are devoted to the matrix elements that correspond to the production 
mechanisms illustrated by the diagrams in Fig. |l|, namely, (a) the one-body term and the 
two-body terms that arise either from (b) pion rescattering or (c) the exchange of heavier 



mesons. In Sec. [ID, we collect a number of calculational details and give the expression for 



the total cross section. 

A. One-body matrix element 

The one-body term [Fig. |l|(a)] can be viewed as a nucleon radiating a pion in the dis- 
torting potential of the other nucleon. The interaction between the two nucleons is needed 
in order to conserve energy and momentum. 

We start with the pseudovector interaction Lagrangian between the pion field (j) and the 



nucleon field ip 13 , 
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Ant = -^^757^r7A-9^0, (1) 

where fi^NN is the pion couphng constant and m-,, is the pion mass. A pseudoscalar La- 
grangian (proportional to ip'y5Tip ■ 0), instead of Eq. (|l]), would yield the same one-body 
contribution; however, pseudovector coupling is favored for reasons that are discussed in 

Sec. irm. 



The field operators in Eq. involve sums over operators that create or destroy particles 
of a given momentum |T^. The term 

Mfioc (7570^0 - 757 " ^Wsip) , (2) 

is responsible for the emission of a neutral pion of momentum g by a proton of momentum p. 
Here f/ is a proton Dirac spinor, and p' = p — q is the final proton momentum. The second 
term in Eq. has the nonrelativistic limit 757 —cr, where cr stands for the conventional 
2x2 Pauli spin matrices. This term, which represents the familiar nonrelativistic cr ■ q pion 
coupling, is odd under spatial reflection of the pion, q —>■ — q. It is thus responsible for the 
production of pions with odd angular momentum (p-wave or higher). Since, near threshold, 
the production of p-wave pions is suppressed by the angular momentum barrier, we are 
interested here only in the first term in Eq. (Q), which, being even under spatial reflection, 
can describe the production of s-wave pions. The existence of this term is a consequence of 
Lorentz invariance, which requires both a timelike and spacelike part to be present in the 
original pseudovector {■j^'j^q^) coupling. 

Alternatively, coupling to an s-wave pion can also be generated from a nonrelativistic 
cr ■ q coupling term by invoking Galilean invariance. Clearly, the coupling should involve a 
dot product between cr and the momentum of the pion relative to the nucleon, rather than 
just the momentum q of the pion. However, since is so much smaller than the mass M 
of the nucleon, the relative momentum differs from q by only a small recoil term involving 
the nucleon momentum p. Thus we expect a cr ■ p coupling term from the nucleon recoil, 
which can produce s-wave pions, that is smaller than the cr ■ q term by a factor of m-^/M. 
The nonrelativistic limit of the first term in Eq. (0), 
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TT ( '\ TT ( \ ^nXl'(^ ■ (P + P')Xs /o\ 

U s' {p HloqoUs (p) ^ — , (3) 

is precisely this recoil contribution, where, near threshold, we take go — "^tt- 

It is straightforward to evaluate the contribution of the one-body term in Fig. |l](a) for 
the partial wave ^Pq ^ -^Sq. The result for the corresponding matrix element is 

^--i^r*^°(T) "»('■' (^^9 

Here Uq and Ui are distorted coordinate-space wave functions for the ^5*0 and ^Pq channels, 
p and p' are the initial and final relative momenta of the two nucleons, and q' is the pion 
momentum in the final pp center-of-mass frame. The distorted waves are normalized such 
that 

UL{r) sin (^pr - ^ + , (5) 

as r — * oo, where Sl is the corresponding phase shift. An analogous prescription applies 
when, asymptotically. Coulomb wave functions are used. 

Equation agrees with the result of Ref. 0, except for the extra factor jQ{q'r/2), which 
represents the s-wave pion wave function. Near threshold, the pion momentum q' is small 
and this factor is close to unity. However, the integral in Eq. involves only scattering 
waves and no bound states. It is thus very long ranged and is found to converge only if the 
pion wave function is included. In the work of Koltun and Reitan, the pion wave function 
was neglected, and therefore it appears that their numerical result for Ji (Table II of Ref. 0) 
is in error. 

When the lUCF Cooler data became available, it was immediately clear that the one- 
body cross section, calculated with Ji alone, does not reproduce the data, underestimating 
the experiment by about a factor of five [||]. This discrepancy is so much larger than any 
uncertainty in the calculation of Ji that the conclusion is unavoidable that there must be 
significant additional contributions to s-wave pion production, over and above the one-body 
term discussed so far. 
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B. Pion rescattering 



An obvious mechanism for pion production that involves both nucleons is the so-called 
pion rescattering diagram shown in Fig. |l](b). We note that in Eq. the wave function for 
a free pion is used, while the two nucleons appear as distorted waves. One might argue that 
the pion rescattering diagram represents the first correction that arises from the distortion 
of the pion. The matrix element J^^ for the rescattering contribution, derived from a simple 
phenomenological pion-nucleon s-wave interaction, is given in Ref. [H as 



J jr 



pp' Jo V 2 / 



dr r 



ui{r) . (6) 



The radial function / is defined as /(r) = e ^^/r, where /i = J jrn 



The pion-rescattering matrix element in Eq. (|^) scales with the parameter Ai. This 
parameter is obtained from the appropriate isospin average of the pion-nucleon s-wave scat- 
tering lengths ai/2 and 03/2, corresponding to isospin ^ and |: 



Ai = — ^(^1/2 + 203/2) 



(7) 



Koltun and Reitan [@], as well as other authors more recently use the older value of 



A™ = 0.005 



(8) 



while Lee and Riska |12|, making use of the computer code SAID [15], favor the value 



AI 



LR 



-0.0023 



(9) 



In both cases, pion-nucleon phase shifts have been extrapolated down to threshold from 
energies where scattering data are available. The two values differ because of new data 
added during the past 25 years and also because of different constraints in the extrapolation 
procedure. On the other hand, there exists experimental information that is more directly 
related to the scattering lengths. For instance, from the branching ratios of the decay of 
pionic hydrogen |TH], one obtains ai/2 — 03/2 = 0.263 ± 0.005m~^; from the measured Is 
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width of pionic hydrogen [0, one obtains 2ai/2 + 03/2 = 0.258 ± 0.012m^^. Combining the 
two results yields 

Ai = 0.001 ±0.002 . (10) 

We conclude that so far there is no experimental evidence that Ai differs from zero or, 
consequently, that pion rescattering contributes to pp —* ppn^. Nevertheless, we will explore 
the sensitivity of the calculated cross section to a possible nonzero value of Ai in Sec. |T|. 

C. Exchange of heavy mesons 

In this subsection we discuss the exchange of mesons heavier than the pion. Clearly, me- 
son exchange, where the intermediate state is a positive-energy nucleon, is already contained 
in the distorted waves used for evaluating the one-body term in Fig. |l](a), as is explicitly 
manifest in the construction of the Bonn potential. However, a virtual, negative-energy 
state as shown in Fig. |l](c) is not contained in the nonrelativistic wave functions. Therefore, 
this contribution must be explicitly added to the pion production operator as a two-body 
correction. 

The matrix element for Fig. |l](c) is calculated using elementary Feynman rules [|14| . Let 



us first consider scalar-isoscalar a-meson exchange (other heavy mesons will be discussed 
below). Neglecting distortions for a moment, this matrix element is given schematically by 



'vv\ vv 



7570^0 + 7570% 



2M / ' ' " ' ' " \ 2M 



Ui , (11) 



where is the momentum of the exchanged a meson, and Ui and U[ {U2 and U2) are the 
initial and final Dirac spinors for the first (second) proton. In Eq. (0) we have only included 
the antinucleon contribution VV /2M for the Feynman propagator of the intermediate nu- 
cleon, taking the nonrelativistic limit kQ + Ek — 2M for the energy denominator (there is an 
imphcit sum over the spin indices of V and V). Also, we have only kept the 757o5'o part of 
the pion-nucleon vertex, since we are interested in s-wave pions. The two terms in Eq. (|TTp 
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correspond to the two possible orderings for the emission of the a meson and the emission 
of the pion. 

It is a simple matter to take the nonrelativistic limit of Eq. (|Il|) by expanding to lowest 
order in 1/M. This yields 

KA ^ i \ ■ (P + PQXl 
■^^^"^-Mi^^eTk^j 2M • 

Here k is the momentum transfered by the a meson and p and p' are the initial and 
final momenta of the first nucleon. Equation ([12|) then represents an additional two-body 
contribution to the effective nonrelativistic operator that describes pion production. This 
term can be directly compared to the one-body term in Eq. which has a very similar 
form, except that the two-body term contains an additional factor of 1/M and a factor 
g^/{y^ + TTi^) from the a-meson propagator. It is important to realize that the one-body 
and (T-meson exchange contributions [Eqs. (^) and ([12|)] have the same sign. Thus the 



two-body contribution will interfere constructively with the one-body term. 

We now express Eq. (^) in coordinate space by taking the Fourier transform, and then 
calculate the matrix element with nonrelativistic distorted waves. The momenta p and p' 
become gradient operators that act on the distorted waves. These give rise to a factor 
(d/dr + 1/r) as in Eq. @). The meson propagator is transformed to a radial function /^(r) 
that contains the mass and the coupling constant g^- This function has the general form 

Ur) = ^- , (13) 

4:71 r 

where x = {a,6,uj, p}. Finally, the contribution of Fig. |l|(c) from the exchange of the a 
meson (x = a) to pion production becomes 

■'.-^r^nJi^)Mrm(^.^l)Mr). (14) 



pp' Jo V 2 / M \dr r ^ 

Equation (|T^) is identical to Eq. (|^) except for the extra factor of fcr{r)/M. In contrast 
to the situation with Ji [Eq. (]^)], the inclusion of the pion wave function jQ{q'r/2) is not 
crucial in this case because of the short range of faij). 
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As mentioned later, we obtain the A^A^ distorted waves from the coordinate-space version 
of the Bonn one-boson-exchange potential. The boson masses and coupling constants needed 
to construct this potential have been fit to A^A^ scattering data and are listed in Table A. 3 



of Ref. |]T8[. For consistency, we use the same parameters in calculating the exchange 

contributions discussed here. We also adopt the technique described in Eq. (A. 28) of Ref. 

to include monopole form factors at all meson-nucleon vertices, according to the prescription 



9x 



gx 



A2 

X 



(15) 



where fc^ is the transferred momentum and A^ is the cutoff mass (also listed in Table A. 3 
of Ref. PI). 

In the following, we also consider the contributions to the diagram in Fig. ||(c) from the 
exchange of mesons other than the a meson, again using the corresponding parameters from 
the Bonn potential. We will find later that these contributions are small compared to J^- 

Let us begin with the scalar-isovector 6 meson. Its contribution Jg has the form of 
Eq. (|T^ with the appropriate mass ms and coupling gs used in Eq. (|13|): 

.2 ^oo / n'r\ f^fr) /d l' 



ml 



(16) 



For pp ppTc^ the isospin factors are the same for the exchange of either an isoscalar or an 
isovector meson. The contribution from the exchange of a vector-isoscalar uj meson can be 
calculated in a similar fashion. It contributes a term 



pp' Jo V 2 / 



M 



— + -] + 
dr r \2M J dr 



ui{r) . 



(17) 



The contribution from the exchange of a vector-isovector p meson has the form of Eq. (|T 
with the exception of a nonvanishing tensor coupling: 

7p(r) 



pp 



M 



.dr r . 



l + C, \ df, 
2M J dr 



uiir 



(18) 



Here Cp = 6.1 is the ratio of the tensor to vector coupling for the p meson |TB[. Note that 
Eq. ([18|) describes only the p-meson contribution to the diagram of Fig. 0(c). In principle, 
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there are other contributions of the p meson to the axial charge that arise from a vrp current 



11|. However, the vrp current has an isospin factor proportional to the dot product between 



the pion field and fi x r2, and therefore does not contribute to the reaction pp —>■ ppir^ 



Again, we modify Eqs. (p^)-(p^) to include form factors following Eq. (A. 28) of Ref. ||18|| . 

In principle, contributions from the exchange of pseudoscalar vr and r] mesons to the 
diagram in Fig. [^(c) are possible. These contributions would be large if pseudoscalar coupling 
(75) were used in Eq. (|l]). However, there are several arguments in favor of pseudovector 
coupling (757^?^) |ll9| , ^ . For example, pseudoscalar coupling implies an unrealistically 
large s-wave pion-nucleon scattering length. In contrast, pseudovector coupling yields zero 
(to lowest order) for the scattering lengths, in agreement with experiment and consistent 
with the (experimentally supported) value of Ai = for the pion rescattering term (see 
Sec. [II B|) . After adopting pseudovector coupling, MEC contributions from vr and r] mesons 



vanish to lowest order in 1/M. 

D. Calculational details and total cross section 

The first step in the calculation is the evaluation of the A^A^ distorted waves Ui and Uq 
in the entrance and exit channels. To this aim we make use of the coordinate-space version 



of the Bonn one-boson-exchange potential, as described in Appendix A. 3 of Ref. |18|. The 



form factor correction is applied according to the prescription in Eq. (A. 28) of Ref. [jTS 
and the Darwin term (proportional to V^) is included following Ref. [^. Boson masses, 
coupling constants, and cutoff masses are listed in Table A. 3 of Ref. for all contributing 
mesons. In that table, two sets of parameters are listed [Bonn potential A (BPA) and Bonn 
potential B (BPB)]; they differ in the contribution of the rj meson and in the value of the 
cutoff parameters. In order to test for the sensitivity of the calculation to the distorting 
potential and to connect to earlier results we have obtained distorted waves also from 
the Reid soft-core (RSC) and Reid hard-core (RHC) potentials [p2|. 



The next step is the evaluation of the one-body contribution Ji [Eq. (^]. This depends 
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only on the A^A^ wave functions and the pion-nucleon couphng constant. For the latter, 
we adopt the value /^atat/^^ = 0.075 which is consistent with the Nijmegen phase-shift 
analysis p3|j2^ . In order to compute the very long-ranged integral in Eq. (^), we integrate 
conventionally from to some rmax and then rotate the contour into the complex plane as 
explained in Ref. This changes an integrand that is oscillating like (sin r)/r into one 



that is exponentially damped. For calculations without the Coulomb interaction, we choose 
'"max to be about 7fm. When the Coulomb interaction is included, r^ax is increased to 100 fm 
in order to be able to use the simplest asymptotic form for Coulomb wave functions in 



the complex plane. Our results are insensitive to the exact choice of rmax- Since the initial 
state has a relatively large momemtum {p ^ 1.9 fm~^), we expect that Coulomb effects in 
the entrance channel are small, even at threshold. Thus, in our calculation, the Coulomb 
interaction only affects the ^Sq final state. 

The next step, in principle, is determining the pion rescattering contribution J^r, as in 
Eq. (^). However, since the rescattering parameter Ai seems to be consistent with zero, this 
step will be omitted except when we investigate the sensitivity of the results to a possible 
deviation of Ai from zero. 

Finally, heavy-meson exchange contributions for a, 6, u, and p mesons are calculated 
using Eqs. (Q) and ([16|)-(|T^). These depend on meson coupling constants and masses and 
form-factor cutoff masses. For these parameters we use the same values that define the 
distorting potentials BPA and BPB, discussed above. For this reason, our calculations with 
BPA and BPB distortions are self consistent. Note that the calculation does not contain 
any parameters that are adjusted to pion-production data. 

The matrix element Jtot for the reaction pp ppn^ is then composed of contributions 
from the one-body term Ji, from pion rescattering (in principle), and from heavy- meson 
exchange currents Jmec^ 

Jtot = Ji + JjT + Jmec ■ (19) 

As described in Ref. [@] , the total cross section is obtained as a phase space integral over the 
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square of a matrix element, 

c^tot = ^f^r^^^ dg'g'VUtotP, (20) 
TTpMml Jo 

where j3 is the lab velocity of the projectile, q' is the pion momentum in the final pp center- 
of-mass system, and p' is the relative momentum of the final protons. 

III. DISCUSSION OF RESULTS AND SENSITIVITY TO INGREDIENTS 
A. Contributions to the matrix element 

We now examine the relative importance of the various contributions to the matrix 
element Jtot [see Eq. (|^)]. To this aim, we evaluate Jtot for the typical values ofp= 1.9 fm~\ 
p' = 0.2 fm^^ and q' = 0.1 fm~^. The results are listed in Table |, which shows that the 
(T-meson contribution J„ is about as large as the one-body term Ji. Furthermore, these two 
contributions are constructive; therefore, the a-meson contribution will increase the cross 
section by about a factor of four! 

The next important contribution is from the u meson (25-35% of the Ji or contri- 
bution). In comparison, all other terms are significantly smaller. In particular, the pion 
rescattering term is small even when the large, older value of Ai = 0.005 is used, and 
the isovector 6- and p-meson contributions are very small and tend to cancel each other. 
For the sake of completeness, we include in Jmec the contributions from all heavy mesons; 
however, it is important to keep in mind that the dominating two-body contribution is due 
to the scalar cr-meson exchange. 

This large contribution [Eq. ([I^) ] can be easily understood in relativistic models. In 
these models [|,^^, the nucleon has an effective mass 

M* = M-S (21) 

from a strong repulsive scalar mean field S. This enhances the lower components in the 
Dirac spinors U. If spinors of mass M* are used, Eq. (P) becomes 
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uAp'; M*)757ogot/.(p; M*) ^ _!!!^^£:^_^±p:)^ . (22) 

Thus a reduction of the nucleon mass enhances the s-wave pion couphng. Let us expand 
1/M* to lowest order in 1/M: 

If one assumes that the scalar field arises from the second nucleon, or 5 = fa{r) [see Eq. (p^)], 
then the second term {S/M), when substituted into Eq. immediately yields the cx-meson 
contribution of Eq. (0). Thus the scalar field of the second nucleon affects the Dirac spinor 
of the first nucleon in such a way that the modified spinor has a larger coupling to an s-wave 
pion. This simple argument also shows that Ji and should add constructively. 

It is nevertheless surprising that the a contribution is almost as large as the one-body 
term, especially since involves an additional factor of 1/M. This is explained by realizing 
that it is not that is large, but Ji that is anomalously small. Near threshold, the two 
nucleons approach each other with a relative momentum p ~ 1.9 fm^^, and yet they must 
almost stop in order to produce a pion. Because it is difficult to mediate such a large 
momentum mismatch through the distortions of Fig. |l](a), Ji is relatively small. On the 
other hand, the a-meson propagator in Fig. ^(c) provides an efficient means of transferring 
momentum. As already pointed out, it is interesting that the cr-meson mass (550 MeV, 
according to Table A. 3 of Ref. |TB[) is comparable to p ^ 1.9 fm~^. 

The interpretation of the cr-meson contribution deserves comment. Clearly the a is 
not a sharp resonance. Instead it is a simple phenomenological model for the important 
intermediate-range attraction in the NN interaction. Assuming that this attraction trans- 
forms as a Lorentz scalar provides a natural explanation of the spin dependence of the NN 
force I^J^. We expect a more complicated model of the intermediate-range attraction (such 
as correlated two-pion exchange) to yield a similar J^- contribution provided the attraction 
transforms as a Lorentz scalar. 
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B. Total cross section results 



In Fig. 1^ our calculation is compared with the available pp —* ppT^^ total cross section 
data as a function of 77, the maximum pion momentum in the overall center-of-mass frame 
in units of (77 = qma.x/fnn), or, alternatively, the projectile energy T in the laboratory. 
The lUCF Cooler data [jT| are shown as solid dots (note that there is a 6.6% uncertainty in 
the normalization of the data that is not shown), while data from previous work are marked 



with crosses squares bars |^, and diamonds |^. As is well known by now, the 
one-body term alone (dashed line) greatly underestimates the data. However, when the two- 
body contributions are included (solid line), the measured cross sections are reproduced to an 
extent that is truly remarkable in view of the fact that none of the parameters of the model 
have been adjusted to pion-production data. The dot-dashed line, which has been obtained 
without the Coulomb interaction, demonstrates that Coulomb repulsion is responsible for a 
fairly sizeable reduction of the cross section near threshold. All calculations shown in Fig. |^ 
use the BPA distorting potential. 

It has been pointed out earlier |^ that the energy dependence of the s-wave cross section 
follows from phase space and the final-state interaction between the two (charged) protons. 
This is sufficient to reproduce the shape of the measured cross section up to r/ ~ 0.6, where 
higher partial waves (which we are neglecting) start to contribute Thus it is only the 
magnitude of the cross section near threshold, represented by a single number, that contains 
nontrivial physics information. Our work shows that heavier meson exchange (mainly of the 
a meson), together with the one-body term, is sufficient to explain the observed magnitude 
of the cross section. Furthermore, there are no alternative explanations at this time. In the 
past, it has been suspected that the role of heavy- meson exchange currents is suppressed by 
A^A^ correlations and form factors at the meson-nucleon vertices. We find that this is not 
the case in the A^A^ system: Because of the simplicity of the present reaction, it is possible 
to explicitly include A^A^ correlations by solving for the full two-body wave function. Even 
with these correlations included, remains large. 
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The dependence of J a on the cutoff mass is illustrated in Fig. ^ (using BPA). Note that 
the form factor in Eq. (|T5|) has a normalization at fc^ = that depends strongly on Ag.: 

9a{^)=9. \, ^ ■ (24) 

a 

Indeed, the primary effect of this form factor is to change the value of 5'o-(0) rather than 
the momentum dependence of the interaction. It is this coupling near A;^ = that plays 
a dominant role when adjusting the parameters of one-boson-exchange potentials to NN 
scattering. Therefore it may be more meaningful to compare results with different A^- values 
at a fixed ga{^) rather than at a fixed g^j. Figure ^ shows only a modest decrease of with 
decreasing A^. (at fixed g^). However, if one keeps 5'cr(0) fixed rather than g„, actually 
increases very slightly with decreasing Ao-. If two potentials fit phase shifts with different 
cutoff masses, we expect approximately similar values of gcr{0) (rather than g^^). If this is 
the case, the Jo- contribution will be almost independent of Ao-. 

We now examine the sensitivity to the various distorting potentials that are mentioned 
in Sec. [11 D[ Figure ^ shows the total cross section, calculated with the Coulomb interaction 
included and without the pion rescattering term, for the two somewhat different one-boson- 
exchange potentials BPA (solid line) and BPB (dotted line), and for the RSC (dot-dashed 
line) and RHC (dashed line) potentials. For the phenomenological Reid potentials there 
is no way to unambiguously determine the meson-exchange contributions. Therefore, we 
simply adopt the meson couplings and cutoff masses from the Bonn potential BPA for both, 
RSC and RHC. This allows us to study the effects of a change in only the distorted waves 
Mo and Ml. The RHC wave function is identically zero at small distances; in contrast, the 
RSC wave function is nonzero, but still small. This enhancement of the wave function at 
small r leads to a modest increase in pion production. At small r, the BPA and BPB wave 
functions are almost identical, but larger still than the RSC wave function, resulting in a 
slightly increased cross section. However, some of the difference between the Bonn and Reid 
results is due to the fact that BPA and BPB generate slightly larger on-shell phase shifts 
as compared to the Reid potentials. This is because the Bonn potentials were fit to pn 
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rather than pp data. Nevertheless, the range in cross section for different potentials is still 
relatively modest. We point out that all calculations without meson exchange (the lower 
four curves in Fig. ^ greatly underestimate the data. 

Based on the small value for Ai in Eq. (p^, we do not expect pion rescattering to 
be important. For completeness sake, we still wish to examine the effect of a possible 
rescattering contribution. All calculations in Fig. ^ are carried out with the BPA and with 
the Coulomb interaction included. The dotted lines correspond to Ai = 0.005, the dot- 
dashed lines to Ai = —0.0023, and the solid line is without rescattering. The effect from 
pion rescattering is relatively small even for the large value Ai = 0.005, and all calculations 
with only one-body and pion-rescattering contributions (the lower two curves in Fig. ^) fall 
substantially below the experimental data. Clearly pion rescattering as described by [see 
Eq. (P)] cannot explain the difference between Ji and the data. However, one must realize 
that Eq. (^ is based on a simple on-shell model for pion-nucleon scattering. It is conceivable 
that the pion-nucleon interaction might be modified due to the fact that the intermediate 
pion in Fig. |^(b) is off the mass shell. This is a topic that deserves further study. 

IV. IMPLICATIONS AND CONCLUSIONS 

We have calculated s-wave pion production in pp —>■ ppn^ by considering the one-body 
term, pion rescattering, and two-body meson-exchange processes. We confirm that the one- 
body term underestimates the data by about a factor of five, and that pion rescattering is 
indeed small. On the other hand, we find a large contribution from the exchange of heavy 
mesons (in particular, the scalar a meson) coupling to the negative-energy state of a nucleon, 
as in Fig. |l](c). This meson-exchange contribution (MEC) is large enough to explain the 
discrepancy between one-body production and the data, and, when taken into account with 
self-consistent distortions, leads to an excellent fit to the data without parameters that are 
adjusted to pion-production information. 

The theoretical description of the pp ppn^ reaction close to threshold is clean and 
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simple. Only a single partial wave is allowed in either entrance or exit channel, pion rescat- 
tering is suppressed, and intermediate A isobars are expected to be unimportant (in marked 
contrast to most other reactions involving pions). Furthermore, the simplicity of the system 
allows a full treatment of A^A^ correlations. 

The operator for the production of s-wave pions has the same form as the axial-charge 
operator. Therefore, we conclude that the axial charge in nuclear systems is much larger 
than one is lead to believe from one-body predictions. This agrees with the conclusion from 



a number of calculations of first-forbidden /? decays in nuclei [|Tll. However, these nuclear 
studies are less conclusive because of structure ambiguities. The pp — > ppir^ system is free 
of such ambiguities. 

The pp — > ppiT^ reaction is sensitive to two-body contributions because of the nature of 
the axial-charge operator and due to a large momentum mismatch that naturally favors the 
exchange of heavy mesons. 

Uncertainties in the meson-nucleon form factors, the two-nucleon wave function, and pion 
rescattering have been examined and are much smaller than the two-body meson-exchange 
contributions. We are thus led to the conclusion that the near-threshold pp ppn^ data 
provide direct experimental evidence for meson- exchange contributions to the axial current. 
This is analogous to the electro-disintegration of the deuteron which provides direct evidence 
for electromagnetic MEG. However, in our system the MEC is a very large effect, compared 
to a ~ 10% contribution to the cross section in the case of the electromagnetic current . 

Furthermore, near-threshold pp — ppvr^ data provide direct evidence for exchange cur- 
rents from heavy mesons. To the best of our knowledge, all previous experimental evidence 
for MEC has involved exclusively pions. This is significant because of the possibility that 
heavy-meson contributions are greatly suppressed by A^A^ correlations and meson-nucleon 
form factors. We have shown that this is not the case in the pp system. The importance of 
heavy-meson exchange currents could be very significant for experiments planned at CEBAF, 
which, involving higher momentum transfer, are likely to be sensitive to such short-distance 
effects. 
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Relativistic nuclear models characteristically feature large Lorentz scalar and vector po- 
tentials [0,^]. In these models, the a meson describes the important intermediate-range 
attraction in the A^A^ interaction and gives rise to a scalar potential that reduces the effec- 
tive mass of a nucleon. This reduction in M* enhances the lower components of Dirac wave 
functions and this increases the axial charge. This change in the Dirac wave functions is also 
what provides the natural description of a large range of nucleon-nucleus scattering data, 
in particular spin observables 0. The fact that we find that this same effect explains the 
observed pion production cross section also provides an indirect experimental confirmation 
of this key feature of relativistic models. 

The Jrj term in our calculation can be viewed more generally as a term involving some 
intermediate-range attraction in the A^A^ interaction that transforms like a Lorentz scalar. 
This contribution does not necessarily have to arise from an elementary narrow a meson. 
Instead, it could well be an effective representation of a more complex mechanism. In any 
case, and whatever the microscopic origin of this Lorentz scalar attraction, its importance in 
the present calculation provides evidence for a large relativistic effect in the A^A^ interaction. 

Future theoretical work on the pp — > ppn^ reaction should be devoted to a study of 
the off-shell aspects of pion rescattering. Also, the present study should be extended to 
include the next higher partial waves, as they become important with increasing bombarding 
energy. This is especially important, since a measurement of the spin-dependent total cross 
section is planned at lUCF that will allow the separation of p-wave pion contributions. This 
provides the data necessary to search for meson-exchange contributions to Gamow- Teller 
matrix elements. Since we have found a large and unexpected meson-exchange contribution 
to the axial charge, we may well speculate about the existence of a similar contribution to 
Gamow- Teller matrix elements. 

This could have important consequences for the pp — >■ de^v reaction and the solar neu- 
trino problem p2[. This reaction is believed to proceed via a ^Si Gamow- Teller 



transition. An enhancement of as little as 15% in the matrix element for pp de^v (from 
an unexpected MEG) would dramatically reduce the disagreement between theory and ex- 
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periment [Q, because, in the standard solar model, the rate of this reaction sensitively 
affects the central temperature of the sun and thus the high-energy neutrino flux. Of course. 



present calculations of tt and p MEG give only a small contribution. Furthermore, MEG 
for this channel are expected to be smaller than the order v/c MEG contributions to the 
axial charge. However, one may still speculate that some unsuspected MEG or other effect 
could be important. Through very accurate pion production data it may be possible to gain 
(indirect) experimental information about the pp de'^v rate. 
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FIGURES 

FIG. 1. Contributions to the reaction pp ppir^ . Shown are (a) the one-body term, (b) the 
two-body term that arises from pion rescattering, and (c) the two-body term that arises from the 
exchange of heavier mesons (x). 

FIG. 2. Total cross section for pp — > ppi:^ as a function of r/, the maximum pion momentum 
in the center-of-mass frame in units of m.,^, or the projectile energy T in the lab frame. The 
recent lUCF Cooler measurement [O is shown by solid dots; data from previous work are marked 



with crosses |27], squares bars and diamonds [30|. The solid line includes heavy-meson 
exchange and the Coulomb interaction, and the dotted lines describe the error bands for the vrA'^ 
coupling constant /^jVAf/'^^ ~ 0.075 ± 0.003 and the rescattering parameter Ai = 0.000 ± 0.002 
(see text). The dashed line neglects heavy- meson exchange and includes the Coulomb interaction, 
while the dot-dashed line includes heavy-meson exchange but neglects the Coulomb interaction. 
(Note that the Coulomb interaction is included in all other curves in this paper.) All calculations 
use BPA. The solid curve is the same as in Figs. ^ and H. 

FIG. 3. Dependence of on the cutoff mass A^- for an initial nucleon momentum of 
p = 1.9 fm~^, a final nucleon momentum of p' = 0.2 fm^^, and a pion momentum of q' = 0.1 fm^^. 
The dashed line is obtained by keeping fixed, while for the solid line 5(j(0) was kept constant 
[see Eq. (p4|)]. The calculations use BPA, which has the standard cutoff Ao- = 2GeV. 

FIG. 4. Sensitivity of the calculated cross section to the distorting potential. Shown is the 
cross section divided by rf' as a function of 77. The data are the same as in Fig. |2|. The curves are 
for the Bonn potentials A (solid lines) and B (dotted line), and for the Reid soft-core (dot-dashed 
line) and Reid hard-core (dashed line) potentials. The solid curve is the same as in Figs. § and 
^. The lower four curves are without Jmec (the curves are labeled in the same manner as above, 
except a dot-dot-dashed curve is used for BPA). 
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FIG. 5. Sensitivity of the calculated cross section to the rescattering contribution using BPA. 
Shown is the cross section divided by rf as a function of r/. The data are the same as in Fig. 
The curves are for Ai = 0.005 (dotted line), Ai = (solid line), and Ai = —0.0023 (dot-dashed 
lines). The solid curve is the same as in Figs. |2| and The lower two curves are without Jmec- 
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TABLES 

TABLE L Contributions to the total matrix element Jtot for an initial nucleon momentum of 
p = 1.9 fm~^, a final nucleon momentum of p' = 0.2 fm~^, and a pion momentum of q' = 0.1 fm^^. 
The Bonn one-boson-exchange potentials (BPA and BPB) and the Reid potentials (RSC and RHC) 
are described in the text. The pion-rescattering contribution J^^ assumes Ai = 0.005 [see Eq. (^]. 



Potential 


Ji 








Js 


Jp 


BPA 


-0.174 


-0.163 


-0.062 


-0.042 


-0.009 


0.009 


BPB 


-0.165 


-0.152 


-0.063 


-0.044 


-0.024 


0.011 


RSC 


-0.143 


-0.135 


-0.050 


-0.038 


-0.006 


0.008 


RHC 


-0.126 


-0.140 


-0.048 


-0.040 


-0.007 


0.009 
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